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Abstract 

qq ■ We study the scaling limits of the L-state Restricted Solid-on-Solid (RSOS) 

t-h , lattice models and their fusion hierarchies in the off-critical regimes. Starting 

with the elliptic functional equations of Klumper and Pearce, we derive the 
' Thermodynamic Bethe Ansatz (TBA) equations of Zamolodchikov. Although 

0^ ■ this systematic approach, in principle, allows TBA equations to be derived 

\ for all the excited states we restrict our attention here to the largest eigen- 

value or groundstate in Regimes III and IV. In Regime III the TBA equations 

are massive while in Regime IV there is massless scattering describing the 

• (i) 
. renormalization group flow between distinct A\ coset conformal field theo- 

j> \ ries. Regimes I and II, pertaining to Zx_i parafermions, will be treated in a 

k^J ' subsequent paper. 



1 Introduction 



Many integrable Quantum Field Theories (QFT) |T[] can be obtained |2| by per- 
turbing Conformal Field Theories (CFT) with respect to an appropriate scaling 
operator. Examples include both massive and massless QFTs. From a microscopic 
point of view, these theories can be obtained as the continuum scaling limit of 
two-dimensional integrable lattice models A primary tool for the study of such 
theories is the Thermodynamic Bethe Ansatz (TBA) (||, [J. For the most part, the 
TBA analysis has been restricted to the groundstate but the method has been 
generalized to treat a few excited states fg, ||. A recent review of these and related 
topics is given by Mussardo |K|. 

^^E-mail address: pap@maths.inu.oz.au 
2 E-mail address: nienhuis@phys.uva.nl 
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In this paper we are interested in a systematic derivation of TBA equations 
starting at the level of the microscopic integrable lattice models. The approach 
developed here should be applicable to all excitations. Although the methods are 
general, we work here in the framework of the Restricted Solid-on-Solid (RSOS) 
lattice models of Andrews, Baxter and Forrester [llj and their fusion hierarchies Jl2| . 



The key to the derivation of the TBA equations is in essence to "solve" the elliptic 
TBA functional equations of Kliimper and Pearce |y§ in the appropriate scaling 



regime for the finite-size corrections to the transfer matrix eigenvalues. The approach 
of the present paper provides an alternative approach to the QFT transfer matrix 



approach introduced recently by Bazhanov, Lukyanov and Zamolodchikov |TJ . 

The layout of the paper is as follows. We begin, in this section, by defining the 
RSOS lattice models and describing their relation to perturbed CFTs in the four 
distinct integrable regimes. We define the commuting transfer matrices and state 
the functional equations of Pearce and Kliimper [p~3 1 . We also introduce the notion 



of the scaling limit and give an overview of the general structure of the TBA equa- 
tions. In Sections 2 and 3 we derive the TBA equations relevant to Regimes III and 
IV respectively. Regimes I and II will be treated in a subsequent paper. The mas- 
sive TBA equations in Regimes II and III have been derived previously by Bazhanov 
and Reshetikhin J7|. Their approach uses the Bethe ansatz equations and proceeds 
subject to a conjecture on the string solutions. Although conjectured by Zamolod- 
chikov H , there appears to be no previous derivation of the massless TBA equations 
in Regime IV. In some sense this is the most interesting regime since it describes the 
massless renormalization group flow between distinct coset conformal field theories. 
We conclude the paper with some discussion of future work. 

1.1 RSOS Models 

The RSOS(p, q) models are restricted solid-on-solid models in which heights on 
the sites i of the square lattice take values in the set {1, 2, 3, ... , L}, with L > 3, 
subject to the additional condition that the values a iy aj of heights on adjacent sites 
of the lattice satisfy the constraints 

< cii < L, < (aj — aj + m)/2 < m, m < ai + aj < 2L — m + 2 (1-1) 

where m = p for a horizontal pair and m = q for a vertical pair. The Boltzmann 
weights W p ' q (u) depend on a spectral parameter u, a nome or temperature variable 
t = p 2 and a crossing parameter A = 7r/(L + 1). These weights are assumed to 
vanish unless the adjacency constraints are satisfied. There are four distinct off- 
critical physical regimes: 

Regime I: -tt/2 + A < u < 0, -1 < t < 

Regime II: -tt/2 + A < u < 0, < t < 1 , , 

Regime III: 0<u<A, 0<t<l ' ' ' 

Regime IV: < u < A, -1 < t < 
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There are two lines of critical points at t — 0, one separating Regimes I and II 
and another separating Regimes III and IV. In particular, the RSOS(jo, p) models 
are conformally invariant on these critical lines. On the Regime III/IV critical line 
the models are described by the coset conformal field theories JU 



with central charges 



3p 
p+2 



1 - 



2(p + 2) 



(L + l)(L + l-p) 



;i.3) 



;i.4) 



On the Regime I/II critical line the models are described by Zl-i parafermion 
theories |16|] with central charges 



2(L-2) 
L + l 



;i.s) 



The face weights of the RSOS(p, q) models, with max[p, q] < L — 1, are con- 
structed by fusing blocks of p x q elementary faces together. The elementary faces 
with (p, q) = (1, 1) are given by the ABF models. In this case the adjacency condi- 
tion reduces to 



1 < cij, dj < L, 



Explicitly, the non-zero face weights W 



u 



— Clj = ±1. 

W(u) are given by 



W 
W 

w 



oil 


a 


a 


a =F 1 


a 


a±l 


a =F 1 


a 


a 


a±l 


oil 


a 



II 



It 



II 



^i(A-m) 



9a±l 



i?i((a- l)A)^((a+ 1)A)\ 1/2 ^(w 



<?a=Fl 

$i(aA ± m) 
t?i(aA) 



tf?(aA) 



*i(A) 



;i.6) 

;i.7) 

;i.s) 
;i.9) 



where $i(u) = $i («,£>) and the g a are arbitrary ^-independent gauge factors that 
cancel out on a periodic lattice and do not effect the eigenvalues of the transfer ma- 
trices. In this paper we use standard elliptic theta functions as given in Gradshteyn 
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and Ryzhik [17] 

0l(lt,p) 



2p 1/4 sin m JJ ( 1 - 2p 2n cos 2m + p 4n ) ( 1 - p : 

71=1 
OO 

2p 1/4 cosw JJ(1 + 2p 2n cos2w + p 4n )(l 

n=l 

oo 

]J(1 + 2p 2n ~ 1 cos 2m + p 2 ^" 1 ))^ - p 2n ) 

n=l 

oo 

]J(1 - 2P 2 "" 1 cos 2m + p 2 ^)^ - p 2n ). 



P^) 



(1.10) 

(1.11) 

(1.12) 
(1.13) 



n=l 



We express the nome p in terms of a real parameter e > by 



P 



Regimes II and III, 
Regimes I and IV, 



(1.14) 



so that t = p 2 = ± exp(— 2tte). In particular, the $1 functions satisfy the quasiperi- 
odicity properties 

^i(m + tt,p) = -A{u,p) (1.15) 
^(u-iln^p) = -p-V 2 S?i («,£). (1.16) 

To define the RSOS(p, q) face weights it is convenient to fix a particular gauge 

g a = (-l) a/2 y/A(aXj. (1.17) 
With this choice of gauge the RSOS(p, q) face weights are given by 

a • (/.• ISA 



,,l Clq+l b q+ l 



ai hi 



u) = flwu)- 1 E II W 

' k=0 ai,...,a„k=l 



if O-k+1 bk+1 
O-k bk 



independent of the values of the edge spins b 2 , . . . ,b q . Here 



p-i 

n 

j=0 



(1.18) 



(1.19) 



and the p x 1 face weights are given in turn by 

/ bk bk+i 



\ai a p+ \ 



u 



<i2,... ,a p k=l 

independent of the edge spins 62, • • • ,b p . 



Ofc dk+l 



u + (k — p)X 



(1.20) 



1.2 Transfer Matrices and Functional Equations 



Consider the RSOS(p, q) models and suppose that a and b are allowed height con- 
figurations for two consecutive rows of an iV column lattice with periodic boundary 
conditions. Then the elements of the row transfer matrices of the RSOS(p, q) models 
are given by 



N 







•) 




a j+1 





;i.2i) 



where a^ + i = a\ and bx+i = b\. The Yang-Baxter equations 



^W p " 



f 9 
a b 



u W p > ; 



e d 
f 9 



e 9 
f a 



v)W p ' s ( 9 I 
a b 



u + v )W q ' s [ d g C b 
a -v) W p ' q 



e d 
9 c 



u 



;i.22) 



imply that, for fixed p, these transfer matrices belong to commuting families 

T p ' q (u)T p ' q \v) = T p ' q \v)T p ' q (u). (1.23) 

Starting with the fusion hierarchy of Bazhanov and Reshetikhin |Tj|, Kliimper 
and Pearce have shown that the commuting transfer matrices satisfy the inversion 
identity hierarchy 



T p ' q T p ' q = f p J p I + T p ' q+1 T p ' q -\ 1 < q < L - 1 



where 



T p k q = T p > q (u + k\), f p = s p (u) N . 



These functional equations close with 



Tg- 1 = 0, T p '° = f p _J, T P ' L = 0. 



If we further define 



rpp,q-lrpp,q+l 
.p,q _ ± X 1 



l<q<L-2 



;i.24) 



;i.25) 



;i.26) 



;i.27) 



then Kliimper and Pearce [0 have also shown that these equations can be recast 
in the form of TBA functional equations 



(I + ^ l )(I + tl 



1 < q < L - 2 



X28) 
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where 

t P,o = t p,L-i = Q (1 2g) 
From quasiperiodicity, it follows that 

— p-^ — = doubly periodic, = doubly periodic (1.30) 

J{q-2)/2 

with period rectangle in the complex u plane given by 

[ ii x me, Regimes II and III 
period rectangle = < (1-31) 
I 7r x 2nie, Regimes I and IV. 

In Regimes I and IV there is an additional symmetry within the period rectangle 

t p ' q (u ± tt/2 + me) = t p ' q (u). (1.32) 

These statements apply to the entries of the transfer matrices and to their eigenval- 
ues. 

At the critical point with t = and in the braid limit 

lim = (1-33) 

U— >±iOO 

the TBA functional equations reduce to a simple recursion relations which can be 
solved immediately to give 

, P , q _ sin(#) sin((g + 2)0) _ sin 2 ((g + 1)0) 

sin sin 



with the quantization 



= -^, m, = l,2,... ,L. (1.35) 



1.3 Scaling Limit and Perturbed CFTs 

Let a denote the lattice spacing. Then the (continuum) scaling limit of the RSOS 
models is given by 



where 



N -> oo, a -> 0, t -> 0, /i = jV|t|" = fixed (1.36) 



4 ' 



l2(L- 1) 
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Regimes III and IV 
, Regimes I and II. 



(1.37) 



Here v is the critical exponent associated with the correlation length of the 
RSOS(p,p) model which is in fact independent of the fusion level p. The expo- 
nent v is related to the usual specific heat exponent a by the hyperscaling relation 
2 — a = dv. Strictly speaking, this relation breaks down in Regime IV for which 
v' = 2u. It is this fact that leads to Regime IV being massless. Also, in some 
cases for p > 1, the free energy is analytic |12] so a is not defined. We ignore these 
problems and scale precisely as indicated in ( [L.36| ) and (|1.37|) . 
In the scaling limit we introduce variables 

At u 

R = lim Na, m = lim — (1.38) 

N—>oo, a^O t^O, a^O a 

where R is the radial distance and m is a mass. It follows immediately from these 
definitions that 



4/i = mR. (1.39) 

For the RSOS(p, p) models in Regimes III and IV, the off-critical elliptic solution 
corresponds to perturbation away from the conformal critical point by the operator 
$ with conformal weights 

A = A = ^^. (1.40) 
L + l K J 



For the RSOS(l, 1) models this is the thermal operator $ = $i 



3- 



1.4 Zamolodchikov's TBA Equations 

The Thermodynamic Bethe Ansatz (TBA) equations describe the scattering of n 
types of relativistic particles on a circle of radius R. Here we will be interested in 
the TBA equations introduced by Zamoldchikov 0. All of these take the form of a 
system of coupled nonlinear integral equations 

n i roc 

e^e) = rmRviie) - ^ — / d9' ^(9 - 6') log(l + e~^) (1.41) 
j=i 

with the Casimir energy 

n /»oo 

E(R) = -Y^ 1 ^ 1 ^z/,(fl)log(l + e-^ (e) ). (1.42) 
j=i 

Here m ; are the particle masses, 9, & are rapidities, e,^) are quasiparticle pseudo- 
energies, Vi{9) are energy functions and the symmetric kernel satisfies &ij{9) = 



7 



<&ji(9). For diagonal scattering, which ensues in Regimes I and II but not Regimes III 
and IV, the kernel is related to the S matrix by 

d_ 
89 



% = log,%(0). (1.43) 



Typically, the energy functions are given by v{6) = |e ±£l or u{6) = cosh^. 

In the scaling region, the finite-size corrections to the largest eigenvalue T(u) are 
related to the Casimir energy by 

sin 1/ I 
-togT(«) = Nf(u) + ^ r E(R) + o(-) (1.44) 

where f(u) is the bulk free energy and the anisotropy angle is given by 

{(L + l)ix, Regimes III and IV 
2{L + l)u (1-45) 
— , Regimes 1 and 11. 
Li — 1 

Despite the appearance oil/N corrections, the system is not in general conformally 
invariant. The system, however, will be conformal at critical points. Such critical 
points can occur in the infrared limit (R — > oo or /x — > oo) in addition to the critical 
point in the ultraviolet limit (R — > or /x — > 0) with 

RE(R) -> -vrc/6 (1.46) 

where c is the central charge. 

The scaling behaviour is manifest in the behaviour of the zeros of T(u) in the 
complex u plane. The imaginary period of T(u) is either 2ne or 4ne where ire = 
2(logiV — log/i)/(L + 1). For large N and finite nonzero /x, the scaling behaviour 
occurs in the two regions 

Im(#) = ±(0 + logiV-log^) (1.47) 

where the rapidity 8 is independent of N. This equation relates the rapidity 9 to 
the spectral parameter u in these two scaling regions. 



2 Regime III 

In this section we derive the TBA equations in Regime III. The strategy is to "solve" 
the TBA functional equations in the scaling limit for the finite-size corrections to 
the eigenvalues t p ' q (u) of the transfer matrices t p,q (u). In these equations the depen- 
dence on the temperature t enters only in the combination t v which we replace with 
fi/N. Regarding /x as fixed there only remains a dependence on N, but since these 
functional equations are exact for finite N they contain all the information required 
to calculate the finite size corrections. To do this we follow closely Kliimper and 
Pearce ]13| . 
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2.1 Integral Equations and Finite-Size Corrections 

As in Kliimper and Pearce we consider the analyticity strips of t p,q (u): 

P-9-2 ( p-q + 2 
A < Re(w) < A. 



(2-1) 



In Regime III the leading contributions to the eigenvalues t p ' q (u) in its analyticity 
strip are given for large N by 



t p ' q (u) = { 



''const 



fP,1 
const i 



; M(L + l)u/2, t v ) 
! t? 2 ((L + l)«/2,t") 



u\-\N 



q = p 



p^q 



(2.2) 



where the constants are given by 

' sin(gcr) sin((g + 2) a) 



sin 2 a 



^const = \ 4 cos a cos r, 



sin((g — p)t) sin((g — p + 2)t) 
sin 2 t 



and 



a 



p + 2' 
m'-TC 

L + l-p' 
For the largest eigenvalue, we have 



m] = 1,2, 



1,2, 



7T 



7T 



1 < g < p — 1 
9 =P 

p + 1 < q < L 

... ,p+l 
... ,L-p. 

7T 



(2-3) 



L + l 



a 



p + 2' 



1 — p 



(2.4) 
(2.5) 

(2.6) 



Other choices will be needed for the excited states. Note that the leading contribu- 
tions to the eigenvalues t p,q (u) give the correct periodicity, have the required zero 
and pole structure and that the constants are such that the functional equations 
(|1.28|) are satisfied including in the scaling limit. 

We introduce functions of a real variable by restricting the eigenvalue functions 
to certain lines in the complex u plane 



A q (x) 



l + a q (x) 



(2.7) 
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Here and the sequel we will usually suppress the dependence on p. Let us also 
introduce finite-size corrections terms £ q (x) by writing 



a q {x) = e q {x)£ q {x) 



with 



e q (x) 



.A{ix/2,t v ) 
''•d 2 {ix/2,V) 



1 N 



, q = p 
q^p 



The TBA functional equations then take the form 

£ q {x - 7Ti/2)£ q {x + m/2) = A q - 1 (x)A q+1 (x). 



(2.8) 



(2.9) 



(2.10) 



The functions £ q (x) and A q (x) are analytic and nonzero in the strip — tt < Im(x) < 7r 
with period (L + 1)tt£. We represent their logarithmic derivatives by Fourier series 



LI 



[\og£ q (x)}' = L\eT^ 

k=—oo 
(L+l)ne/2 



1 



(L + 1)7T£ J-IL+i 



[log^(x)]'e dx 



(2.11) 
(2.12) 



(L+1)tt£/2 



with similar equations relating the Fourier coefficients A\ to A q [x). 

So taking the logarithmic derivative of ( |2.10| ) and equating Fourier coefficients 
gives 



e (£+i) £ _|_ e (i+x)e 



and hence 



[bgf (*)r= 



e( L + 1 ) £ 



k=— oo 



g(L+l)e _|_ g (£+1)6 



^r 1 + 



4- 



(2.13) 



(2.14) 



Substituting the integral expression for A q k (x) and evaluating the sum on k gives the 
nonlinear integral equation 



[\og£ q (x)]' 



(L+l)7re/2 



dy k(x - y) {[\ogA q -\y)}' + [\ogA q+ \y)}'} . (2.15) 



'-(i+l)Tre/2 

After integration this yields 

log a 9 = log e q + k * log A^ 1 + k*log A q+l + D q 



(2.16) 
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where * denotes convolution (over the given interval) and D q are integration con- 
stants. 

The kernel is given in terms of standard Jacobian elliptic functions and theta 
functions by 



fc=_oo e( i+1 ) £ + e ( L + 1 > £ 



ILa ( 2I ' x ~\ 1 gUMOM^g r 2 i7i 

#2(0, ^3(0,^3(^,4') 



27r-i9 2 (ia;, g') 

where / and J' respectively denote the complete elliptic integrals of nomes 

..2 

q = e -TLTWe } q> = e -*(L+i> = ±L. (2.18) 

In the critical limit, when e — > 00 and <f — > 0, the kernel simplifies to 

lim k(x) = k(x) = — (2.19) 

e->oc K J v ' 2vrcoshx V ' 



To evaluate the constants D q we set 



1 



(L+1)7T£ 



{a) = TiTWsl a{x)dx (2 ' 20) 



Then integrating (|2.16|) and using the result 



I 



to 

we find that 



(L + 1)tT£ 1 

k(x)dx = - (2.21) 



(log^) = l - ((A q " 1 ) + (A 11 - 1 )) + D q . (2.22) 

Comparing with (|2.10 ) we see that, for each q, D q is a multiple of iri related to the 
branches of the logarithms and winding. In particular, for the largest eigenvalue, 
there is no winding and D q = 0. Alternatively, we could work directly with the 
Fourier series of the logarithms rather than the logarithmic derivatives. In this case 
we would need to allow for multiples of 2ni on the right side of ( [2.13 ) with the same 
end result. 

Let us factor the eigenvalues T p,q {u) into bulk and finite-size correction terms 

T p > q (u) = T^(u)T^ te (u) (2.23) 



11 



then as in Kliimper and Pearce 

%(«)%(« + A) = /(« - A)/(« + gA) (2.24) 

and 

T^ te (n)T^ te (u + A) = 1 + (2.25) 

where now T™ ltc {u) is a doubly periodic function of u. It is again useful to introduce 
functions of a real variable by restricting the eigenvalue function to certain lines 

= T^ tc [j^x + ^|^A) (2.26) 

so this relation becomes 

b q {x - m/2)b 9 (x + m/2) = A q {x). (2.27) 

Since b q and A q are non-zero and periodic we can solve as before by introducing 
the Fourier series of the logarithmic derivatives to obtain 

log 6" = k *\og A q + C q (2.28) 

where C q are constants and the convolution is over the finite interval. Explicitly, 
the finite-size 1/N corrections ( 2.28 ) for q = p are 

r(L+l)we/2 

logW(x)= dyk{x-y)log{l + a p {y)) + C p (2.29) 

J-(L+l)ne/2 

We evaluate the constant as before by integrating ( 2.29Q to obtain 

(&P) = + C p . (2.30) 

Comparing with ( |2.27| ) we see that, for each p, C p is a multiple of ni so it does not 
contribute to the 1/N corrections. 



2.2 TBA Equations and Casimir Energy 

So far our integral equations are completely general and are exact for finite N. 
They could be used, for example, to study the exponentially small corrections to 
the eigenvalues off-criticality. In this subsection, however, we specialize the form of 
these equations appropriate to the scaling limit. 

The system size N only enters the integral equations through the function e p (x). 
This function is periodic in x with period 2 (log N — log /i). For iV large the function 
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is exponentially small in N except in the two scaling regions when x is of the order 
of log N or — log N. Let us set 



Z 2 = e -( x + l °g N ) 



N 



r 



N' 



(2.31) 



Then in these scaling regions we find 



loge±(x) 



lim loge p (±(x + log N)) 

N— >oo 



lim log 

N— >oo 



1 + t 2v z 2 



1 + t lv z 



1 + Z 2 

= -2{e~ x + fi 2 e x ). 
We assume that the functions a q and A q scale similarly and set 



a±(x) 



Al(x) 



lim a p (±(x + logN)) 

N— >oo 

lim A p (±(x + logiV) = 1 + a q ±(x). 

N—*oo 



(2.32) 

(2.33) 
(2.34) 



Scaling x and y in the same way, the integral equations ( |2.16|) take the following 
simplified form in the scaling limit 



log a q = \oge q + k* log A' 1 ' 1 + k * log A q+1 + D q 
where we suppress the subscripts ± and 

0, q j^p 



log e q (x) 



-2(e x + /i 2 e x ), q — p. 



(2.35) 



(2.36) 



Clearly, this reduces back to the case of Klumper and Pearce when /i — > 0. Indeed, 
the modified equation is obtained just by replacing — 2e~ x by — 2{e~ x + fi 2 e x ) in 



If we now set 



6 = x + log fi 



(2.37) 



and 



a q ( x ) = a q {6- log /i) = e"^ (e) (2.38) 
then — 2(e~ x + fj, 2 e x ) = —4ficosh8 and we obtain the TBA equations (|1.41|) 



00 ^ logq + e-^')) + j_ r ^iog(i + e-^ e 'y) 



cosh(0-0') 2n 
+ D q = mR cosh 68. 



cosh(fl - 6') 



(2.39) 
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where q — 1, 2, . . . , L — 2 and we identify 

4/i = mR. (2.40) 

The Casimir energy is given by the scaling limit of the finite size correction. The 
integral in ( |2.28| ) is over one period (L + l)ne = 2(logiV — log/i). For N large, the 
integrand is of the order o(l/N) unless y is in one of the scaling regions where y is 
of the order of log N or — log N. Hence 

^ rlog TV-log fi 

loglf(x) = - dyk(x-y-\ogN)log(l + a p (y + logN)) 

2 J- (log TV-log/.) 



log TV— log /i 



+ - / dyk(x + y + \ogN)log(l + a p (-y-\ogN)) + C p 

2 J -(log TV-log fj.) 
1 f°° 

dy (e x ~y + fx 2 e~ x ^) log(l + d p + (y)) (2.41) 

+ r dy {e ' x ~ y + ^ ex+v) iog(i + h ~ {y)) + ° p+o o) 



Here we have used the fact that 



k(±(x + logiV)) = e ~^f eX +o[jj). (2.42) 



For the largest eigenvalue we have a + (y) = (i-(y) and 

l°gl£L(«) = ^/_^2/(^ + ^V)log(l + a p (y)) 



oo 

oo 



2 fx cosh x 
Rsin(L + l)u 



N 

where the Casimir energy is 



/oo 
c^coshfllog(l + e- ej,(e) ) (2.43) 
-oo 

E P (R) 



EJR) = / ^cosh^log(l + e-^ (9) ) (2.44) 

2tt y_ 00 

and we have used 4 fx = mR. In the isotropic case u = 7r/2(L+l) and sin(L+l)w = 1. 



3 Regime IV 

In this section we derive the TBA equations for Regime IV. We use precisely the 
same scaling as for Regime III with t < 0. It turns out that in this regime, with one 
exception when p = (L — l)/2, the resulting TBA equations are massless. 
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3.1 Integral Equations and Finite-Size Corrections 

In Regime IV the analyticity strips are the same as for Regime III. The leading bulk 
contributions to the eigenvalues t p,q {u) in its analyticity strip are given for large N 
by 



+p,g 

''bulk 



(«) 



''const 



const 



: M( L + 1 >/2,\t\ 2v 
l ^2({L + l)u/2,\t\ 2 " 



N 



M(L + l)u/2 + (L-2p- 1)tt/4, \t\ 2v ) 
& 3 ((L + l)u/2 + (L — 2p — 1)tt/4, \t\ 2 ») 



-i N 



V 



q = L — p — 1 



const ) 



otherwise 



(3.1) 



where the constants are the same as in Regime III. Again this has the required 
periodicity and zero and pole structure. This applies for p < (L — l)/2. For the case 
p > (L — l)/2, the zeros and poles of order iV can be obtained from the previous 
case using the relation 



+ p,q 



4.P ,g 
V+i' 



p — L — 1 — p. 
In the marginal case p = p' = (L — l)/2 we find 

N 



(3.2) 



fP,q 

''const 



M(L + l)u/2,\tn 



q = p 



# 2 ((L + l)u/2,\t\») 

otherwise. 

As in Regime III we introduce finite-size correction terms £ q (x) by writing 



(3.3) 



fP,g 

k const ) 



a q {x) = t p ' q 
Here, for p < (L - l)/2, 



L + l 



x + 



-A 



e 9 (x)^(x) 



(3.4) 



e 9 (a;) 



A(ix/2, 


t 


2v\ 




t 


2v^ 



'^{ix/2, \t\ 2v ) 



N 



N 



V 



q = L — p — 1 
otherwise 



(3.5) 
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and, for p — (L — l)/2, 

j .1?i(lX/~Z, ^ s 

e q (x) 





t 


n 




t 





' q = P (3.6) 
otherwise 



Following the same Fourier series analysis as in Regime III and allowing for the 
change in periodicity leads to the integral equations 

log a q = log e q + k * log A"' 1 + k * log A q+1 + D q (3.7) 

In this regime the kernel is given by 



j x 

HX) = 2(L + 1W E 



2/kx 



6 2(L+1)£ 



2(L + 1)7T£ fe f-^ g2(£+l) S _|_ g 2(L+l)e 

^(0,^3(0,^3(^,2') 



27r^ 2 (ix,g') 

with the elliptic nomes 



(3.8) 



g = e"W*, g' = e - 2 ^ L+1 ) £ = ^. (3.9) 
In the critical limit, the kernel again simplifies to 

lim k(x) = k(x) = l —— . (3.10) 

e^oc w w 27rcosha; v 7 

The constants can be evaluated as in Regime III, allowing for the different period, 
by setting 

Again D q is a multiple of ni with .D 9 = for the largest eigenvalue. 

Repeating the analysis of the finite-size corrections for Regime IV leads to the 
result 



/(L+l)ire 
dy k{x-y)\og{l + a p {y)) + C p (3.12) 
■(L+1)7T£ 



where C p is a multiple of 7ri that does not contribute to the 1/N corrections. The 
only difference with Regime III is in the period. 
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3.2 TBA Equations and Casimir Energy 

In the scaling regions we find 

loge± = lim loge p (±(x + log N)) 

N— »oo 

-2e~ x , p<(L-l)/2, q = p 

-2/i 2 e x , p<(L-l)/2, g=L-l-p (3.13) 

-2(e" x + /iV), q = p=(L-l)/2. 

Passing to the scaling limit of the integral equations ( |3.7| ) thus yields 

loga g = loge 9 + k * logi^ 1 + k * logA q+1 + D q (3.14) 

where again we suppress the subscripts ±. Next, introducing rapidity variables as 
before we obtain the TBA equations. For p = (L — l)/2 we obtain the same massive 
TBA equation as in Regime III. For p < (L — l)/2, however, we obtain the massless 
TBA equations (|L4lD 



2tt > /_ 00 cosh(0-0') 27t7_ 0O cosh(0-0') 

+ D q = X -mRe- Q S pq + ^mtfe* 5z_i_ Pl9 (3.15) 

where g = 1, 2, . . . ,L — 2 and we again identify 4/x = mi?. 

Turning to the finite-size corrections we see that the integral in ( |3 .12 ) is over one 
period given by 2{L + l)7re = 4log(N/p,). Hence 



/iog(Jv//i) 
A;(x — y — log — ) log(l + a p (y + log — )) 
-log(jV/^) I 1 H 

dy fc(x + y + log — ) log(l + a p (-y - log — )) + C p 

■ log(JV/At) A* A 1 

/oo 
dye*-"log(l + a^(y-log^)) (3.16) 
-oo 



JL 

ttN 



+ r dy e ~ x ~ vl °^ 1 + *-(y - lo ^)) + CP + (a?) 



Here we have used the fact that 



fc ( ± (, + log^)) = ^ +0 (i|. (3.17) 
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2/x cosh a; 
Rsm{L + l)u 

w 



For the largest eigenvalue we have a + (y) = a_(y) and 

2 cosh r f°° 

logT^Ju) = J Jy e-y log(l + a*{y)) 

/oo 
dOe' 6 'log(l + e- £ " w ) (3.18) 
-oo 

# P (i2) 

Hence the Casimir energy is 

EJR) = -— / d6e- e \og(l + e'^ e) ) 

= - — / ^r e - 9 log(l + e-^ (9) ) + e e log(l + e-^- 1 ^ (e) )l (3.19) 

where we have identified 4 fx = mR and used the symmetry between e p {6) and 
€L-i-p{— 0). In the marginal case when p — {L — l)/2 we see that e p (^) is even and 

777 /"°° 

E P (R) = / ^cosh^log(l + e-^ (9) ) (3.20) 

in agreement with the massive case in Regime III. 

4 Discussion 

In this paper we have given a systematic derivation, at least for the largest eigen- 
value, of the TBA equations for the RSOS lattice models and their fusion hierarchies 
in the off-critical Regimes III and IV related to A± coset models. Interestingly, in 
the case of Regime IV, this appears to give the first derivation of Zamoldchikov's 
massless TBA equations describing the renormalization group flow between distinct 
coset theories. In a subsequent paper we will similarly, derive the massive TBA 
equations for Regimes I and II pertaining to Zu-x parafermions. In this case the 
TBA equations for the largest eigenvalue in Regimes I and II turn out to be exactly 
the same due to a duality between the two regimes. 

The systematic derivation of TBA equations introduced in this paper seems to 
afford many advantages over the approach based on Bethe ansatz equations. First, 
and most importantly, the analysis in this paper should in principle generalize to 
allow for the treatment of all the excited states. Second, the present methods can also 
be extended to treat systems with a boundary. We hope to explore these possibilities 
fully in future work. Already, for the case of in the massless Regime IV with 
fixed boundaries, it is possible |L9| to obtain the TBA equations for all excited states 



and to give a complete classification of these eigenvalues. This enables a complete 
mapping of the flow of eigenvalues from the tricritical Ising to the critical Ising 
conformal fixed points. 
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